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The total entropy production generated by the dynamics of an externally driven systems exchang- 
ing energy and matter with multiple reservoirs and described by a master equation is expressed as 
the sum of three contributions, each corresponding to a distinct mechanism for bringing the system 
out of equilibrium: nonequilibrium initial conditions, external driving, and breaking of detailed 
balance. We derive three integral fluctuation theorems (FTs) for these contributions and show that 
they lead to the following universal inequality: an arbitrary nonequilibrium transformation always 
produces a change in the total entropy production greater or equal than the one produced if the 
transformation is done very slowly (adiabatically). Previously derived fluctuation theorems can be 
recovered as special cases. We show how these FTs can be experimentally tested by performing 
the counting statistics of the electrons crossing a single level quantum dot coupled to two reservoirs 
with externally varying chemical potentials. The entropy probability distributions are simulated for 
driving protocols ranging from the adiabatic to the sudden switching limit. 

PACS numbers: 



INTRODUCTION 



In statistical mechanics, thermodynamic laws are 
recovered at the level of ensemble averages. The past 
decade has brought new insights into nonequilibrium sta- 
tistical mechanics due to the discovery of various types 
of fluctuation relations valid arbitrarily far from equilib- 

rium i, i, i, i, i, 0, B i, EO, im, E m, Q [H, [ii, [13 . 

These relations identify, at the level of the single re- 
alization of a statistical ensemble, the "trajectory 
entropy" which upon ensemble averaging reproduce the 
thermodynamic entropy. They therefore quantify the 
statistical significance of nonthermodynamic behaviors 
which can become significant in small systems [l8|, [l9| . 
Various experimental verifications of these FTs have 
been reported [2(J, [2l|, [H, [H, [H, [H, . 
In this paper, we consider an open system, described 
by a master equation (ME), exchanging matter and 
energy with multiple reservoirs. The system can be 
externally driven by varying its energies or the different 
temperature or chemical potentials of the reservoirs. 
There are three mechanisms for bringing such a system 
out of equilibrium: preparing it in a nonequilibrium 
state, externally driving it, or putting it in contact with 
multiple reservoirs at different temperatures or chemical 
potentials thus breaking the detailed balance condition 
(DBC). We show that each of these mechanism makes 
a distinct contribution to the total entropy production 
(EP) generated by the nonequilibrium dynamics of the 
system. The two first contributions are nonzero only if 
the system is not in its steady state and are therefore 
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called nonadiabatic. The third contribution is equal 
to the EP for slow transformation during which the 
system remains in a steady-state and is therefore called 
adiabatic. We derive three FTs, for the total EP and 
its nonadiabatic and adiabatic contribution and show 
that they lead to exact inequalities valid arbitrary far 
from equilibrium. Previously derived FTs are recovered 
by considering specific types of nonequilibrium trans- 
formations. Steady state FTs [l2|, |l^ are obtained 
for systems maintained in a nonequilibrium steady state 
(NESS) between reservoirs with different thermody- 
namicproperties. The Jarzynski or Crooks type FTs 
0, 0, are derived for systems initially at equilibrium 
with a single reservoir which are externally driven out 
of equilibrium by an external force. The Hatano-Sasa 
FT [To!, [Hi] is recovered for externally driven systems 
initially in a NESS with multiple reservoirs. 
To calculate the statistical properties of the various 
contributions to the total EP and to demonstrate the 
FTs, we extend the generating function (GF) method 
[6] to driven open systems. Apart from providing clear 
proofs of the various FTs, this method is useful for 
simulations because it does not require to explicitly 
generate the stochastic trajectories. Some additional 
insight is provided by using an alternative derivation of 
the FTs similar to the Crooks derivation j^, Q , where the 
total EP and its nonadiabatic part can be identified in 
terms of forward-backward trajectory probabilities. By 
doing so, we connect the trajectory approach previously 
used for driven closed systems [8|, 0, H with the GF 
approach used for steady state systems p|, [l2| . 
We propose to experimentally test these new FTs in 
a driven single orbital quantum dot where the various 
entropy probability distributions can be measured by 
the full electron counting statistics which keeps track of 
the four possible types of electron transfer (in and out 
of the dot through either lead). Such measurements of 
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the bidirectional counting statistics have become feasible 
recently We calculate the entropy probability 

distributions, analyze their behavior as the driving is 
varied between the sudden and the adiabatic limits, and 
verify the validity of the FTs. 

In section |TT1 we present our stochastic model and in 
section IIIII we describe the various contributions to the 
total EP generated during a nonequilibrium transforma- 
tion and the inequalities that these contributions satisfy. 
In section |lVl we define the various trajectory entropies 
which upon ensemble averaging give the various contri- 
butions to the EP. We then present the GF formalism 
used to calculate the statistical properties of these tra- 
jectory entropies. In section [Vl we derive the various FTs 
and the implied inequalities. Alternative derivations of 
FTs in terms of forward-backward trajectories are given 
in appendix [Aj By considering specific nonequilibrium 
transformations, we recover most of the previously de- 
rived FTs. Finally in section [VH we apply our results to 
the full counting statistics of electrons in a driven quan- 
tum dot. Conclusions are drawn in section [Vlll 



II. THE MASTER EQUATION 



each satisfying 
W^"} (A) 

,(A) ' 



exp <^ pj,{\)[{emW - e^/(A)) 



(4) 



If all reservoirs have the same thermodynamic prop- 
erties (temperature and chemical potential /i), the 
steady state distribution coincide with the equilibrium 
distribution p^(A) = p^{X) which satisfies the detailed 
balanced condition (DBC) 

<L (a)pl^, (A) = w!::i^ix)p^^{x) . (5) 

As a consequence of dH and ([5]), the equilibrium distri- 
bution then assumes the grand canonical form 

e,..^ ^ exp[-/?(A)(e^(A)-MA)iV„»)] 

-Pm V^J "(A) ' 

where S(A) is the grand canonical partition function. 
However, in the general case where the reservoirs have 
different /^^^^ and /i^^^ , the DBC does not hold and p^(A) 
is a NESS. 



We consider an externally driven open system exchang- 
ing particles and energy with multiple reservoirs. Each 
state m of the system has a given energy and Nm par- 
ticles. The total number of states m is finite and equal 
to M. The probability to find the system in a state m at 
time t is denoted by Pm{t). The evolution of this proba- 
bility is described by the ME 



m' 

where the rate matrix satisfies 

^H^„^,m'(At)=0. 



(1) 



(2) 



We assume that if a transition from m to m' can occur, 
the reversed transition from m' to m can also occur. 
Various parameters, such as the energies Cm of the 
system or the chemical potential jij, and the temperature 
(3~^ of the V reservoir can be varied in time externally 
according to a known protocol. This is described by the 
dependence of the rate matrix on several time-dependent 
parameters A^. If the transition rates are kept constant, 
the system will eventually reach the unique steady state 
solution p^(A) which satisfies p^(A) = [1^. 

The transition rates will be expressed as sums of con- 
tributions from different reservoirs v 



Wrr 



(3) 



III. THE ENTROPIES 

The Gihhs entropy of the system is a state function 
defined as 



S{t) = - ^Pm{t) Inp^(t) . 

m 

Using ^ and ([2]) , the system EP reads 
S{t) = -^Pm{t)\np^{t) 

m 

= - E Wm,m'{><t)Pm'{t)\n 



(7) 



Pm{t) 



(8) 



' Pm'{t) 

This can be partitioned as [H, [H, [sS, [H, [H 

S(t) ^ Stot(t) - Sr(t) ^ (9) 
with the total EP 

Stotit) ^ - Yl w^ir,L(AtW(t) (10) 



In — > 

W^il'(AtW(i) 

and the reservoir EP (also called medium entropy or en- 
tropy fiow) 

sAt)^- y: <L(A.wwin ^ Tjf ; . (11) 



W""' A\t) 
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We note that Stot{t) > follows from wj^^^,{Xt)pm'{t) > 
if m' 7^ m and Inx < x — 1 for x > (if 
m' = m the log in zero), by using the fact that 
E«,.Vr(r,;«(A,)p^(A,) = and E«,.W^ir,L(At) = 0. 
S{t) is the contribution to Stot{t) coming from the 
changes in the system probability distribution and 
Sr{t) is the contribution coming from matter and energy 
exchange processes between the system and its reservoirs. 

We further separate the reservoir EP into two compo- 
nents Ofii 



Sr{t) = Sea:{t) ^ Sa{t) , 

with the excess EP 



(12) 



(13) 



Sex{t) = - XI ^m,nf (At)Pm'(*)ln- ^ 

= Xp«(i)lnp^(At) 



and the adiabatic EP (also called housekeeping entropy 

[a El) 



Sa{t) 



(14) 



m.ra' ,u 



m > . 



<1'(A0p^'(a*) 



The positivity of ([Ej) follows from the same reason as 
([To]). If a transformation is done very slowly, the system 
remains at all times in the steady state distribution 
Pm{t) = p^(At). Such a transformation is called adia- 
batic. We then have Stot{t) = Sa{t) and Sex{t) = —S{t). 
We also notice that Sa{t) = when the DBC is satisfied. 



We next define the state function quantity 

Pm{t) 



Sb{t) = - Y,Pm{t) In 



(15) 



which is obviously zero when the system is at steady 
state. When considering a transformation between 
steady states, AS'5(T,0) = drSbir) = Sb{T)-Sb{0) = 
0. We call Sb{t) the boundary EP (the terminology will 
be explain shortly) and separate it in two parts 



Sb{t) = Sna{t) - Sd{t) , 

where the nonadiahatic EP is 

Pm{t) 



Sna{t) = -'^Pm{t)ln 



p'A{>^t 



(16) 



(17) 



and the driving EP is 

Sd{t) = Ep™(*)'^™(a*) 



(18) 



with 



■lnp™(At 



(19) 

The positivity of (pTj) is again shown in the same way 
as for ([To]) and ([T4|) . If no external driving acts on the 
system, A is time independent and from ([T8|) . Sd{t) = 0. 
For an adiabatic transformations, since Pm{t) = p^(A^), 
from ([T5|) and ([TTj) . we see that Sna{t) = as well as 
Sb{t) = 0. From ([T6|) . this also means that Sd{t) = 0. 
Therefore, Sd{t) 7^ only for nonadiabatic driving. Using 
(fT7|) with ([8]) and ([13]), we find 

Sna{t) = Sex{t) + S(t) = Stot(t) - Sa(t) . (20) 

It is clear from the last equality why we call Sna{t) the 
nonadiabatic EP. The inequality Sex{t) > —S{t) which 
follows from the first line is a generalization of the "sec- 
ond law of steady state thermodynamics" [HI, [sl, [s^ 
derived for transitions between steady states. 
We next summarize our results 

Stot{t) = Sna{t)^Sa{t)>0 (21) 
Sna{t) = Sd{t)^St{t)>0 (22) 
Sa{t) > 0. (23) 

The total EP is always positive and can be separated into 
two positive contributions, adiabatic (which are nonzero 
only when the DBC is violated) and nonadiabatic effects. 
The latter can be due to a nonadiabatic external driving 
acting on the system or to the fact that one considers 
transformation during which the system is initially or 
finally not in a steady state. We therefore have a mini- 
mum EP principle stating that the total EP for arbitrary 
nonequilibrium transformations takes its minimal value 
if the transformation is done adiabatically (very slowly). 
The equality sign in ([2T]) is satisfied for adiabatic trans- 
formations which occur at equilibrium. The equality sign 
in ([22]) holds for adiabatic transformations. The equality 
sign in ([23]) only occurs when the DBC is satisfied. 



IV. TRAJECTORY ENTROPIES 

The evolution described by the ME can be represented 
by an ensemble of stochastic trajectories involving 
instantaneous jumps between states. This will allow us 
to define fiuctuating (trajectory) entropies. 



A. Definitions 

We denote a trajectory taken by the system between 
t = and t = T by m(^) = {0 — mo ^ mi ^ • • • rnj-i 
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rrij • • • rriN-i ttin — T}. At t = the system is 
in mo, and stays there until it jumps at ri to mi, etc., 
jumps at Tat from rriN-i to tun and stays in mAr until 
t = T [see Fig. [1]. A/" is the total number of jumps during 
this trajectory. 
We next introduce various type of 'trajectory entropy 



m2 / / mj-2 



mj / / niN- 



FIG. 1: Representation of a trajectory m(^). 

production' (TEP). We will see at the end of this section 
that when ensemble averaged, these correspond to the 
various EP defined in section Hill 
The trajectory Gibbs entropy is defined as 



(24) 



where Pm^^^it) represents the value of Pm{t) along the 
trajectory m(^). The system TEP is given by 



s\m 



Pm{t) 



Pm{t) 



N 



Pn 



PmArj) 



(.25) 



The first term represents the smooth changes of s[m(^), t] 
along the horizontal segments of the trajectory on Fig. 
[T] during which the system is in a well defined state. 
These changes are only due to the time dependence of 
the probability to be on a given state. The notation \m^^) 
means that the m in the expression changes depending on 
which horizontal segment along the trajectory one con- 
siders. The second term represents the discrete changes 
of s[m(^), t] along the vertical segments of the trajectory. 
These changes are singular and only due to the change 
in the system state. 

Separating the trajectory system TEP similarly as the 
system EP in section HIH we get 



Stot[m(^r),t] =<5[m(^),t] +5r[m(^),t] 
where the total TEP is 

. r n _ Pm{t) 

Stot[m^r),t\ 



Pm{t) 



N 



+ y5{t-Tj)\n 



(26) 
(27) 



AK,)Pm,{Tj) 



and the reservoir TEP 

N 



yy-mi ,771-,- 



w. 



(28) 



We further separate the reservoir TEP into 

Sr [m(^) ,t]= Sa [m(^) , t] + Sex [^(r) , t] , (29) 

with the adiabatic TEP 

N 

S [ , '] 

j=l " '"^j-ii'if^j \-' j /J-rrij \-' j 

and the excess TEP 



4m,^„t] ^Y.'it- r,) m (30) 

VVmj-i,mj y^Tj jPrrin V^r.- ) 



Sex[m^r).t] = l^S{t - Tj)\n-- . 

j=l Prnj-A^Tj) 

The nonadiabatic TEP 

Pm{t) 



Sna[mi^r)^t] 



Pmit) 



(31) 



(32) 



+ }_^5{t- Tj) In , . 

is made of the sum of two terms 

Sna [^(r) ,A=^d [^(r) , ^] + ^fo [m(^) , t] (33) 

the boundary TEP 



Pm{t) 



^</>m(,)(At) 

dXt 



Pm{t) 

Pm,-,(r,)p^,.(Ar,) 



(34) 



and the driving TEP 



5rf[m(^),t] = Xt- 



P%_Mrj)Pm,{rj) ' 

{Xt) 



dXt 



As in section Uni since 

Sna[m^r)iA = s[m(^) , t] + 5ea^ [m(^) , t] , (35) 

we get 

Stot [m^r) , A = ^na [^(r) , ^] + [m(^) , t] . (36) 

We generically denote these TEP by a[m(^),t]. The 
change of a[m(^),t] along a trajectory m(^) of length T 
is given by 



Aa[m(^),T] 



dta[m(^^^^t\ . (37) 



Notice that /S.s[m(^^^^T] and A55[m(;^),T] are state func- 
tion TEP 

A.[m(.),T] = In = s^AT) - SmM , (38) 

where Sm{t) = — Inp^(t), and 

A55[m(,),T] = {Sm^{T)-^m^{XT)) (39) 

-(Smo(O) - ^mo(Ao)) • 

The other TEP are path functions. 



5 



B. Deriving trajectory entropies from measured 
currents 

Using dl]), the reservoir TEP can be expressed as 

Sr[m^r).t] = (40) 

where the heat current between the u reservoir and the 
system is 

N 

and the matter current between the u reservoir and the 
system 

AT 

The currents are positive if the system energy (matter) 
increases. 5j,^jj.{t — tj) is a Dirac distribution centered 
at time tj only if the transition is due to the reservoir 
V = Vi. Otherwise it is zero. We thus confirm that the 
reservoir EP is the entropy associated to system-reservoir 
exchange processes. 

We assume that the parametric time dependence of the 
energies, temperatures and chemical potentials is known. 
Except in degenerate cases for which two different transi- 
tions between states have the same energy difference and 
number of particle difference, the trajectory of the sys- 
tem can be uniquely determined by measuring the heat 
and matter currents between the system and the reser- 
voirs. The system steady state probability distribution 
can be calculated by recording the steady state currents 
for sufficiently long times for different values of the en- 
ergies, temperatures or chemical potentials. The driven 
system probability distribution can in principle be cal- 
culated by reproducing the measurement of the currents 
multiple times. All trajectory entropies containing the 
logarithm of the transition rates can be expressed in term 
of a combination of the reservoir EP (directly measurable 
via current) and other trajectory entropies which can be 
expressed in term of the system probability distribution 
(actual or steady state). Therefore, provided the cur- 
rent measurements can be repeated often enough to get 
a good statistics, all the trajectory entropies are in prin- 
ciple measurable. 



C. Statistical properties using generating functions 

The GF formalism allows to compute the probability 
distributions and all statistical properties of the TEP 
without having to generate the trajectories themselves. 



It further provides a direct means for proving the FTs. 
The GF associated with the changes of a[m(^), t] along a 
trajectory is given by 



G(7,t) = (exp{7Aa[m(^),t]}) , 



(43) 



where (•) denotes an average over all possible trajecto- 
ries. The probability that the system follows a trajectory 
with the constraint A A = Aa[m(^),t] at time can be 
obtained from the GF using 

P(AA,t) = ((5(AA- Aa[m(^),t]) (44) 
By inverting (|44|). we get 

/CO 
dAAe''"^'^P{AA,t) . (45) 
-oo 

The moments of the distribution are given by derivatives 
of the GF 

/A fcr .n d''G{'y,t) 



1,2, 



• (46) 



7=0 



In order to compute the GF, we recast it in the form 
G(7,i) = E5-(7,i), (47) 



where 



(48) 



is the product of the probability to find the system in 
state m at time t multiplied by the expectation value 
of exp {7Aa[m(^), t]} conditional on the system being in 
state m at time t. Since Aa[m(^),t] = for a trajectory 
of length t = 0, we have Qmin-, 0) = Pm(0). We also have 
G(0,t) = land gm{^,t)=p^{t). 
The time derivative of (|47j) gives 



(49) 



where gmij^t) depends on the TEP of interest. Below, 
we will derive equations of motion for the ^rn(7,^)'s as- 
sociated to the various TEP. 



1. State function trajectory entropy production 

The generating function associated with a state func- 
tion TEP Aa[m(^),t] = am{t) — an(0) like A5[m(^),t] 
or Asb[m(^^^^t] may be straightforwardly obtained using 
dSj) with gHI). We get 

Gil^t) = ^exp{7(a^(t) - an{0))}pm{t)Pn{0) . (50) 
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2. Excess trajectory entropy production 



5. Total trajectory entropy production 



-^ecc [^(r) 7 ^] acquires an amount 5ea:(m,m') = 
In (p^(A^)/p^/(At)) each time a transition from a state 
m' to m occurs and it remains constant along a given 
state m of the system. This means that 



(51) 



+ exp {7Sex(m, m')}(exp {^/\Sex[m(r),t]})m' , 
which using (|38|) can be rewritten 



'(A05i'f^(7,t).(52) 



Each time a transition along the system tra- 
jectory occurs, Stot[n^{T)^t] acquires an amount 

^^{W^l'i^t)Pm'{t)/W^ljXt)pm{t))- In addition 
it also changes by an amount —d(^lnpjn{t))/dt during 
an infinitesimally small time on a given state m of the 
system. Combining the two, we have 



■Jtot) 



E 



Pmjt) 
Pm{t) 



(tot) 



(56) 



W^lMPmit) 



3. Reservoir trajectory entropy production and currents 

Each time a transition along the system trajectory 
occurs, 5r[m(^),t] acquires an amount Sr{m^m') = 

^^{W^U^{Xt)/W^)^^{Xt)). Similarly to the excess en- 
tropy, we get 

^H^,(A,)x7 



(7,*) = E 



w'-i (At) 



<l'(At)5L1(7,t) (53) 



This equation has been used in the study of steady state 
FTs 0, [13, . 



In (|4Q|) . we expressed the reservoir TEP in terms of 
currents. The time integrated individual currents give 
the heat and matter transfer between the v reservoir 
and the system (?ieL[^(r), ^] = /o ^^'4eL[^(r), ^1 and 
^mat[^(r),^] = dt' I^^l^[m^^),t']. Their statistics can 
be calculated using 



'(At))} 



exp {7^1. (^m - A^rnO}^i>(A.)^rn'(7,^) ,(^4) 

where 7 is a vector who's elements are the different JheatS 

and 7^at's. The GF calculated from ([54|) is therefore as- 
sociated with the join probability distribution for having 
a certain heat and matter transfer with each reservoir. 



4. Adiabatic trajectory entropy production 



6. N on- adiabatic trajectory entropy production 

Like the total TEP, Sna[^{T)->A acquires an amount 
In (p^(At)p^/(t)/p^,(At)pm(^)) each time a transition 
from a states m' to m occurs, and also changes by an 
amount —d(^hiprn(t))/dt during an infinitesimally small 
time on a given state m. This gives 



flL""H7,t) = -7 



Pm{t) 



~.{na) 



E 



Pm{t) 
f^{Xt)Pn.'{t) -^-' 
Pt'{^t)Pm{t) 



(57) 



7. Driving trajectory entropy production 

Since A5^i[m(^),t] exclusively accumulates along the 
segments of the system trajectory, we get 



■Jd) 



7<Pm{Xt)gi^H7.t) (58) 



It follows from (|46|) that the average change of a TEP 
is obtained from its GF by differentiation with respect to 
7 at 7 = 0. By differentiating the GF evolution equations 
of this section one recover the evolution equation for the 
EPs of section IIIII The EPs are therefore the ensemble 
average of the TEPs introduced in this section A{t) = 
(d[m(;^),t]) and AA{T^O) = {Aa[m(^^^^T]). 



Each time a transition along the system trajectory 
occurs, acquires an amount Sa{m^m') = 

ln(H^ir,L'(At)p^'(At)/H^ir,^,^(At)p^(At)). We therefore 
get 



E 



(55) 

^'^^^Vir,L(A.)e)(.,.). 



m' ,m\ 



V. FLUCTUATION THEOREMS 
A. General integral fluctuation theorems 

We can easily verify that gmij = —^^t) = Pm{t) is 
the solution of the evolution equations ([56|) and ([57j) . It 
immediately follows from probability conservation and 
dm) that Gto^(-l,^) = Gna(-1,^) = 0. Summing both 
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side of over m, we also verify that Ga{—l,t) = 0. 

( z) 

Because 1,0) = Pm(0), Gz{—1^0) = 1 where z = 

tot^na^a. Therefore, we find that Gz{—l,t) = 1. Using 
(|43|) . this results in the three FTs 

(exp{-A5tot[m(^),t]}) = 1 , (59) 

{exp{-Asna[m^r)it]}) = 1 ^ (^0) 
(exp{-As,[m(,),t]}) = l. (61) 

These FTs hold irrespective of the initial condition 
and the type of driving. Using Jensen's inequality 
(e^) > e^^\ they imply the inequalities (f2T]) - ([23l) . 

Eq. ([59|) is the generalization to open systems of the 
integral FT for the TEP obtained earlier for closed sys- 
tems [l4]. The TEP 1^ needs to specify which reservoir 
is responsible for the transitions occurring along the 
trajectory (by labeling the rates with reservoir index). 
This point, made earlier for open system at steady state 
is generalized here for driven systems with an arbi- 
trary initial condition. Eq. (|6Q|) will be shown in next 
section to reduce to the integral Hatano-Sasa FT [IT] for 
systems initially in a steady state. Eq. (|6T]) generalizes 
the integral FT for the adiabatic entropy [l5| previ- 
ously derived for closed system initially in a steady state. 

Some additional insights can be gained by an alterna- 
tive proof of the FTs ([59|) and (|6Q|) which use a forward- 
backward trajectory picture of the dynamics. This is 
given for completeness in appendix [Al 



B. Transitions between steady states 

We consider a system initially (t=0) at steady state 
and subjected to an external driving force between t = tdi 
and t = tdf' For tdi > t > 0^ the system remains in the 
steady state corresponding to A = At^.. The time pro- 
tocol of A^ during the driving tdf > t > tdi is arbitrary. 
If ttr is the characteristic transient time needed for the 
system to reach a steady state from an arbitrary distri- 
bution, for t > tdf -\- ttr, the system is in the new steady 
state corresponding to A^^^ . The system is measured be- 
tween t = and t — T. 



1. Fluctuation theorem for the reservoir entropy production 

We restrict our analysis to cases where the system is 
at steady state at t = and the driving starts at least a 
time ttr after the measurement started: tdi > Ur- 
We use the braket notation where \p(t)) is the probability 
vector with components Prn{i) and Wt denotes the rate 
matrix. |/) denotes a vector with all components equal 
to one. The ME ([1]) now reads 



The generating function for the reservoir EP ([53|) for 7 = 
— 1 evolves according to the adjoint equation of (|62|) 

\g^^\-l,t))=wl\g^^\-l,t)). (63) 

The initial condition of (|62]) and (|63]) is |p(0)) = 
|^/(^)(-l,0)) = \p^\\u,)). The formal solution of ^ 
for t < tdi before the driving starts reads 

Gri-l,t) = {I\exp{VVhW^) (64) 

We now insert a closure relation in term of right and 
left eigenvectors of the adjoint rate matrix between the 
evolution operator and the initial condition. Because the 
rate matrix and its adjoint have the same eigenvalues (all 
negative and one zero), for ttr < t < tdi, only the right 
and left eigenvector associated with the zero eigenvalue 
survive. Since the right [left] eigenvector of is |/) 
[{Pm{^td^)\]^ we get for ttr <t< tdi 

Gr{-l,t) = {I\e^p{m}\I){f'\f') . (65) 

For longer times, even when the system starts to be 
driven, |/) remains invariant under the time evolution 
operator as can be seen using ([2]) in (|63|) . We get 

Gr{-hT) = M{p''\p'') for T> ttr, (66) 

where M = (^|^) is the total number of states. This 
implies the following integral FT for the reservoir TEP 

M 

M > (exp{-As.[m(,),T]}) = M ^ {p^iXtjf > 1(67) 

The equality on the l.h.s (r.h.s) is satisfied if p^(At^.) = 
^n,m (Pmi^tdi) = 1/M). Jeuscu's inequality implies 
ASr{T,0) > 0. Note that since As[m(^),T] is a state 
function, it is easily verified that 

M 

M J2 ^'-(O) = (exp{As[m(,),T]}) . (68) 

m=l 



2. The Hatano-Sasa fluctuation theorem 

We assume that the driving starts at the same time 
or later as the measurement {tdi > 0). We define T' > 

tdf + ttr- 

We have pointed out at the end of section IIIII that for 
transitions between steady states AS^iT' , 0) = 0, so that 

ASna{T\0) = ASd{T\0) = ASd{tdf,td^) > (69) 

We used the fact that ASd{T'^0) starts (stops) evolving 
at tdi {tdf)- The same is true at the trajectory level, since 
from (|39|) we have A55[m(^),T'] = and therefore 



\pit)) = ^t\pit)) 



(62) 



As„a[m^r),T'] = Asd[m(,),T'] , (70) 
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where 

A..K.,r'] . Em^j^ (71) 

AT 

i=o 

The integrand in the third hne contributes only during 
the time intervals between jumps provided the driving 
is changing. Therefore if without loss of generality we 
choose the measurement time such that T > tdf (if T < 
tdf one can redefine tdf as equal to T), Asd[m(^^^^T] = 
As dim (^^^^T']. This means that for a transition between 
steady states, the FT ([6Q|) reduces to the Hatano-Sasa 
FT \XU 

{exp{-Sd[m^r).T]}) = l. (72) 

Alternatively, ([72]) can be proved from ([58]) by show- 
ing that when 7 = -1, g^\-l,t) exp {-(/)^(A^)} = 
p^(At) is solution of ([58|) . 

The FT ([72]) holds for an arbitrary driving proto- 
col. Let us consider the two extremes. For an adia- 
batic (infinitely slow) driving, the inequality in ([69]) be- 
comes an equality. In the other extreme of a sudden 
driving, where = Xt^^ + 9(t - tdi){Xtdf - ^tdJ and 
At = S{tdi){Xt^f - At^J, 

Asd [m(^) , T] = (j)mo (Kf ) - ^mo (K. ) (73) 

becomes a state function and its average takes the simple 
form 

ASd{T,0) = Y.p':,{\u.){M\u,) - 4>^{K.)) ■ (T4) 

m 

Using ([35]) with ([70]) . and since 

As[m(^),T] = (t)mT{Kf) - ^mo(AtaJ , (75) 
we find that 

^Sex [m(^) , T] = 4>rno {M^f ) " </>mT {^t^f ) (76) 

also becomes a state function and its average becomes 

mo ,mT 

(^mo(Atd/) — ^mT (At.,)) .(77) 



C. Transitions between equilibrium states 

For a system coupled to a single reservoir (or mul- 
tiple reservoirs with identical thermodynamical proper- 
ties), the DBC (|5|) is satisfied. A non-driven system in 
an arbitrary state will reach after some transient time 
ttr the equilibrium grand canonical distribution (|6]). We 
again choose T' > tdf + ttr and T >tdf. From the TEP 
of section IIV[ we find in this case 

Asa[m^r).T'] = (78) 
Asr[m(^),T'] = Asex[m^r)iT'] 
Astot[m^r)iT'] = Asna[m^r)iT'] . 

The two FT (|59|) and jBO]) become identical and the FT 
([6T]) becomes trivial. Using ((6|), we also find that Eq. 
([T9]) becomes 

^m{X) = m{^m{X) - fi{X)Nm - 1^(A)} , (79) 

where Q{X) = — InS(A) is the thermodynamic grand 
canonical potential. 

We next consider transitions between equilibrium 
states, so that the procedure is the same as in IV B 21 
but with the DBC ([5|) now satisfied. We therefore have 
Asna[m{T)^T'] = AS(i[m(^), T']. The driving implies 
externally modulating the system energies, the chemical 
potential or the temperature of the reservoir. 

When driving the system energy, using ([71]) and ([79]) . 
we find 

Asd[m^r).T] = Asd[m^r).T'] = Pw[m^,^] - PAn (80) 

where the work is given by w[m(^^^] = J^^^ dtem^^^{Xt) 
and An = n[e{tdf)] - n[e{tdi)]. Both FT, $9ij and ([60]) . 
lead to the same Jarzynski relation [4] 

(exp {-pw[m^r^]}) = exp {-pAQ} . (81) 

When driving the reservoir chemical potential, Eq. 
([TT]) and ([79]) give 

Asd[m^r).T] = Asd[m^r).T']=pw[m^r)]-P^^ (82) 

where w[m(^^^] = — J^^^ dtfiNm^^^ and AQ = Q[jj.{tdf)] — 
n[j^{tdi)]. Both FT, fy) and now lead to 

(exp {-Pw[m^r)]}) = exp {-pAQ} . (83) 

The case where reservoir temperature is driven can be 
calculated similarly. 

D. No driving: steady state fluctuation theorem 

In a NESS, the relations of section ITVl give 

A5na[^(r),^] = ASrf [m(^) , t] = (84) 

As[m(^),t] = -Asex[m^r)it] 

Astot[m^r)it] = A5a[m(^),t] . 
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Furthermore, since AS{t, 0) = 0, we get 
AStotit,0)=ASrit,0)>0. 



(85) 



We shall rewrite the GF evolution equation for the reser- 
voir TEP ([53|) in the bracket notation 



so that 



\g<^^Hj,t)) = Vij)\g^^\j,t)) 



G,(7,t) = (/|exp{V(7)t}|p«*), 



(86) 



(87) 



where |/) is a vector with all elements equal to one. Since 
from ([53|) the generator has the property V(7) = V^(— 7— 
1), its eigenvalues have the symmetry 5^(7) = 7 — 1). 
Furthermore, since exp {V(7)t| is a positive matrix, the 
Frobenious-Perron theorem [291 , Isol . [31] ensures that all 
eigenvalues are negative or zero and that the left and the 
right eigenvectors, |^o(7)) and \£,o{7))^ associated with 
the largest eigenvalue 5^0(7) exist. Adopting the nor- 
malization (^0(7)1^0(7)) = 1, we find for long times 

Grh,t) * exp{sMt}mh))MlW) ■ (88) 
and that 

Gr{-j-l,t) *=°° exp{s5„(7)t} (89) 
(/|eo(-7-l))(eo(-7-l)b^*). 

This means that the cumulant generating function 



1 



lim - In Grijj t) 

t^oo t 



satisfies the symmetry 

F,(7) = F,(-7 - 1) . 



(90) 



(91) 



Using the theory of large fiuctuations this symmetry im- 
plies the detailed steady state FT [1, 



F{-ASr) 



(92) 



where F{ASr) is the probability for a trajectory of the 
system to produce a reservoir TEP equal to ASr- 
Asr[m(^^^^t] grows in average with time because it de- 
pends on the number of jumps along the trajectory. How- 
ever, A5[m(^),t] is bounded. The FT ([92|) can there- 
fore be viewed as a consequence of the detailed FT for 
AStot derived in (|A13p . The long time limit is needed 
in order to neglect the contribution from As[m(^),t] to 

The FT (1571) remains valid at steady state. FTs for cur- 



rents can also be derived 



d at steady t 



VI. ENTROPY FLUCTUATIONS FOR 
ELECTRON TRANSPORT TROUGH A SINGLE 
LEVEL QUANTUM DOT 



We have seen in section HVB I that the various entropies 
can be calculated by measuring the different currents 
between the system and the reservoirs. The counting 
statistics of electrons through quantum dots has recently 

raised considerable theoretical Mjli, EO, Eil, Ea] as weii 

as experimental [23, SB, 46| interest. The single 



electrons entering and exiting a quantum dot connected 
to two leads can be measured. One can therefore 
calculate all currents, deduce the system trajectories 
and calculate the various trajectory entropies presented 
earlier. 

We will analyze the probability distribution for the var- 
ious trajectory entropies in a single level quantum dot of 
energy e connected to two leads with different chemical 
potentials /izy(t), where v = l^r. We neglect spin so that 
the dot can either be empty or filled 1. The ME is of 
the form Q [40|, [13, S, m EO, Ell, El 



where 



Po{t) 



Vt 



Wt 



-Vt Wt 
Vt -Wt 



Vi{t) 
Po{t) 



(93) 



V V 



(94) 



The coefficients characterize the coupling between 
the dot and the lead v with Fermi distribution = 
l/(exp{/3{e-/i,(t)}} + 1). If n = 1, [a] = [energy] = 
[time By renormalizing energies by e, all parameters 
of our model become dimensionless. The steady state 
distribution of the system is 



(st) 



(st) 

Po 



Vt 



Vt ^Wt Vt^ Wt 

and the steady state currents are given by 



{its 



Wt 



Vt + Wt 



Vt 



Vt + Wt 



(95) 



(96) 



We switch the chemical potential of the left lead 
Mo + V{t) using the protocol 



V{t) 



Vf-V^ 



tanh (c(t — tm) + 1) 



(97) 



while holding the right lead chemical potential fixed 
jiir{t) = /io- We can therefore calculate all the trajectory 
entropies' probability distributions using the GF method 
described in section IIV CI We solve numerically the 
evolution equations for the ^^(i7,t)'s associated with 
the different entropies for different values of 7 with the 
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initial condition 5f^(i7,0) = Pm(0). After calculating 
the G(i7,t)'s using (|T7|) , the probability distribution is 
obtained by a numerical inverse Fourier transform (|44j) . 
In all calculations we used /3 = 5,e = l, a/=0.2 and 
ar = 0.1. 

We start by analyzing the different contributions to 
the EP as defined in section Hill for the fii{t) protocol 
shown in Fig. ([2]) a. 

The system is initially in a nonequilibrium distribution 
different from the steady state. The solution of the ME 
([93|) as well as its steady state solution are displayed in 
Fig. [2)3. Between t = and t = 20, is essentially 

constant and the system undergoes an exponential 
relaxation to the steady state. Between t = 20 and 
t = 50, changes from jhq -\- Vi to jiio -\- Vf fast enough 
for the system distribution to start differing again from 
the instantaneous steady state distribution (adiabatic 
solution). After t = 50, fii{t) remains constant and the 
system again undergoes a transient relaxation to the 
new steady state corresponding to /io + Vf. Fig. [Jt 
shows the time dependent EP Stot and its adiabatic 
Sa and nonadiabatic contribution Sna- As predicted, 
these three quantities are always positive [see ([2T]) - ([23|) ]. 
We also demonstrate that Sna only contributes when 
nonadiabatic effects are significant i.e. when the actual 
probability distribution is different from the steady state 
one [pm{t) 7^ Prn\fJ^{t))]. Sa = ouly once at t ^ 33, 
when /j.i{t) = /ir(^) = 0.5 and the DBC is satisfied. 
Otherwise the DBC is broken and Sa > 0. In Fig. [2]i, 
we present the two contributions to the nonadiabatic 
EP Sna^ the driving EP Sd and the boundary EP 
Sb [see (|22|)]. The driving EP Sd only contributes 
when fii{t) changes in time. One can also guess that 
ASb = J^^ dtSb = due to the fact that the change of 
boundary EP during an interval between two steady 
state is zero. ASb = J^^ dtSb ^ because the system 
is initially not in a steady state. Fig. [2^ shows the 
alternative partitioning of the nonadiabatic EP into the 
system EP S and the excess EP Sex [see (|2Q|) ]. Finally 
the splitting of the total EP in the reservoir EP and the 
system EP [see (j9j)] is shown in Fig. [2f. We see that at 
steady state 5* = so that Stot = Sr- 

We next study the statistical properties of the dif- 
ferent TEP for transitions between steady states. The 
probability distributions are obtained using the GF 
method presented in section IIV CI The five driving 
protocols used to change from /io + Vi to /io + Vf 

are presented in Fig. [3l They range from sudden switch 
in (i) to slow (adiabatic) switch in (v). The system is 
always initially in the steady state corresponding to 
Mo + ^ • We will consider measurements which end when 
the system reaches its new steady state at /io + Vf. 
Different measurement times are represented by a,b,c,d. 

In Fig. [31 we display P{AS). Since As is a state 



function, P{AS) is the same for the various protocol. 
Because we consider a two level system. As can only 
take four possible values which correspond to the four 
possible change in the system state between its initial 
and final condition. The transitions 0^0 and 0^1 are 
much more probable because the probability to initially 
find the system in the empty state is much higher 
(0.96) than finding it in the filled state 1 (0.04). The 
transition ^ is more probable than 0^1 because 
the system has a final probability 0.64 to be in its empty 
state and 0.36 to be in its filled state. 

In the left column of Fig. O we depict P{ASd)' Here 
Asb = so that Asd = Asna [see (p!6|) ]. All curves 
(i)-(v) satisfy the FT For the sudden switch (i), 

Asd becomes a state function which only depends on the 
initial state of the system [see ([73]) ]. Asd can therefore 
take two possible values corresponding to the empty or 
filled 1 orbital with a respective probability 0.96 or 0.04. 
When the driving speed slows down in (ii) the peaks are 
broadened. In the adiabatic limit (V), P{ASd) becomes 
a broad distribution with zero average. In the right 
column of Fig. [5l we depict P{ASex)' For sudden switch 

(i) , Asex turn to a state function which only depends 
on the final steady state distribution [see ([76]) ]. It is 
clear from ([76|) that the transitions 0^0 and 1^1 
leads to Asex = and 1^0 and ^ 1 to the same 
Asex with opposite sign. The probabilities to observe 
these transitions follow from the fact that the system 
is initially more likely to be in (prob 0.96) than in 1. 
The probability for the final state (1) is 0.64 (0.36). 
Therefore the most likely transition is ^ followed 
from 0^1. As the driving speed slows down like in (ii), 
the peaks get broadened. Since Asex = Asd — As [see 
([20]) ] and since in the adiabatic switch limit (v) P{ASd) 
is centered around zero, P{ASex) (v) has the same peak 
structure as P{—AS) [see Fig. 0], but broadened by Asd. 

In Fig. [6l we display P{ASr) for different measure- 
ment times and protocols. Plots with same driving but 
different measurement times [(i)a and (ii)d or (ii)a and 

(ii) d or (iii)b and (iii)d] show the evolution of P{ASr) 
in the final steady state. The plots (ii)-(v) satisfy the 
FT (|67|) which for our parameters imply (exp {—Asr}) = 
1.846. The FT is not satisfied for (i) because the driving 
starts at the same time as the measurement [see section 
IVB Ij . To understand the structure of P(A6'r), we time 
integrate (|^Q|) and use the fact that in our model the heat 
current is proportional to the matter current between the 
u reservoir and the system /"^^^ = I^l^ = /^g^^/e where 

N 

&\t) = - T,){Nm, - N^^_,) . (98) 

i=i 

We get 

Asr{t) = -f3j2 f dT(e-fi,{T)\&\T) . (99) 
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In the sudden switch limit (i), we get 



,{T)W''\t). (100) 



Asr{t) = -PY1 



where A/"^^^ = /q drl^^^ (r) is the net number of electron 
transferred from the reservoir v to the system between 
and t. This explains why As^ in (i) only take discrete 
value which are multiples of each other. The distance 
between the peaks of 2.5 observed in (i) is due to the 
right lead only /3(e — /i^) = 2.5 because our parameters 
are such that /3{e — iJiiiT)) = 0]. The new peaks which 
appear in (ii) with a spacing 0.125 are due to the fact 
that the driving starts some time after the measurement 
so that /3{e — /i^(0)) = 3.75 also contributes. As the 
driving speed slows down (iii)-(v), the discrete structure 
broadens and Asr can take continuous values. 

In Fig. [71 we display P{ASa)' All curves satisfy the 
FT (|6T]) . The verification (not shown) is best done on 
the GF {Ga{—l^t) = 1) because the numerical accuracy 
of the tail of the distribution is not sufficient. The 
peak structure of P{ASa) can be understood from 
P{ASr) and P{ASex) because Asr = Asa + Asex- 
This is particularly clear for the sudden switch (i) 
where the possible values of the entropies are strongly 
restricted. Indeed, in (i) each peak of P{ASr) is split 
in three smaller peaks which have the same structure as 
P{ASex)' As the speed of the driving decreases (ii)-(v), 
the peak structure disappears. 

In Fig. El we show P{AStot)' All curves satisfy the FT 
(|59|) [verification was done on the GF (not shown)]. The 
structure of P{AStot) can be understood using P{ASr) 
and P{AS) because Astot = + As. This is clear for 
the sudden switch (i) where the peaks of P{ASr) are 
split in smaller peaks which have the structure of P{AS). 



VII. CONCLUSIONS 

For a driven open system in contact with multiple 
reservoirs and described by a master equation, we have 
proposed a partitioning of the trajectory entropy produc- 
tion into two parts. One contributes when the system is 
not in its steady state and contains two contributions due 
to the external driving and the deviation from steady 
state in the initial and final probability distribution of 
the system. The second part comes from breaking of 
the detailed balance condition by the multiple reservoirs 
and becomes equal to the total entropy production when 
the system remains in its steady state all throughout the 
nonequilibrium process. Both parts as well as the total 
entropy production satisfy an general integral fiuctua- 
tion theorem which imposes positivity on their ensemble 
average. This partitioning also provides a simple way to 
identify which part of the entropy production contributes 



during a specific type of nonequilibrium process [see Fig. 
[9]. Previously derived integral fiuctuation theorems can 
be recovered from our three general fiuctuation theorems 
and in addition we derived a new integral fiuctuation 
theorem for the part of the entropy production due to 
exchange processes between the system and its reservoirs 
(reservoir entropy production). Our results strictly apply 
to systems described by a master equation ([1]). However, 
as has often be the case for previous fiuctuation relations, 
one could expect similar results to hold for other types 
of dynamics. For electron transport trough a single level 
quantum dot between two reservoir with time dependent 
chemical potentials, we have simulated and analyzed in 
detail the probability distributions of the various trajec- 
tory entropies and showed how they can be measured in 
electron counting statistics experiments. 
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APPENDIX A: FLUCTUATION THEOREMS IN 
TERMS OF FORWARD BACKWARD 
TRAJECTORY PROBABILITIES 

We show that the FT ([59|) and (|6Q|) have an interesting 
interpretation in term of the ratio of the probability of 
a forward dynamics generating a given trajectory and 
the probability of the time-reversed trajectory during 
some backward dynamics. This is an alternative to the 
GF approach which connects the detailed form to the 
integral form of the FTs. 

The forward dynamics is described by the ME ([1]) . We 
introduce the probability (in trajectory space) V[m(^^^] 
that the system follows a trajectory m(^) 



dr'W, 



N 



n (Ai) 



exp 



rrij — i ,mj 
T 



i(Ar0j^i7k-i(^-.) 
exp( / drWmN,mN{^r')) . 



The Wm-^rrij-iiKj) factors in this expression represent 
the probability that the system undergoes a given tran- 
sition whereas the exponentials describe the probability 
for the system to remain in a given state between 
two successive jumps. Summation over all possible 
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trajectories will be denoted by 



It consists of We can partition (|A6p in the form 



time-ordered integrations over the N time variables rj 
from to T (this gives the probability of having a path 
with N transitions) and then summing over all possible 
from to oo. Normalization in trajectory-space 
implies that T.^^^^ 7^[^(r)] = 1- 

The backward dynamics is described on the time in- 
terval t = [0, T] by the ME 



(A2) 



where the new rate matrix Wm,m'{Xt) satisfies 
^^Wm,m'{Xt) — 0. We require that the parametric 
time dependence (via the driving protocol At) of the rate 
matrix in Eq. \k2\ is time-reversed compared that of Eq. 
([1]) and that the diagonal part of the rate matrix in Eq. 
(d]) and (|A2]) is the same 



^m,m(At) Wm,m(At) . 



(A3) 



This still leaves room for different choices of Wm,mW\ 
We will later specify two choice of Wm,m'(AT-t) [ (|Alip 
and ()A14p ] that will result in two FTs. 
We define the time-reversed trajectory of m(^) by m(^) = 

• • - mi "^-^^ mo — T}. The probability V[fri{T)\ that the 
system described by ()A2p follows the time-reversed tra- 
jectory m(^) is given by 



V[m^r)\ 



exp 



- PmAr(O) 

T-TN-j + l 

T-rN-j+2 



N 

n 

dr'W, 



(A4) 



mN-j+i,rnN-j+i i^T-r') 

^^XT-TN-j + l) 



niN-j N 



exp 



dT'W„ 



T-Ti 



,(At-x')) , 



where tat+i = T. Normalization in the reverse path 
ensemble implies X]m( ^ ^[^(r)] = 1- 



We consider the ratio of the two probabilities (|Aip and 



In 



(A5) 



Due to ()A3p . the contributions from the exponentials 
(which represent the probabilities to remain on a given 
state) in (|A5p cancel, so that 



r[m(^)] In 



Pmo(Q) 
PmAr(O) 



AT 



In ^ / , 



-.(A.,) 



,(Ar,) 



(A6) 



(A7) 



r 1 1 Pmo(O) , Pm,(Ar,) 

r mr^J = m — ^^^4 + > m — — ^— - — - 

We assume for the moment that r[m(^)] can be expressed 
exclusively in terms of quantities of the dynamics ([T]) i.e. 
a recipe has to be provided to express the tilde quantities 

in (|A6]) [pmN (0) and Wm-li,mj {Kj )] in terms of non-tilde 
quantities. 

In analogy with (|A5p . we define 



f[m(^)] In 



(A8) 



for the tilde dynamics. The previous recipe also 
implies that r[fh(^r)] can be exclusively expressed in 
terms of quantities of the tilde dynamics ()A2p . Eq. 
()A5p together with ()A8p implies that r[m(^)] = — ^[m(^)]. 

The probability F{R) to observe a trajectory such that 
r[m(^)] = during the forward dynamics is related to 
the probability F{—R) to observe a trajectory such that 
f[m(^)] = —R during the backward dynamics 

P(i?) ^ ^P[m(.)],5(i?-r[m(.)]) (A9) 
= E^N(r)]e'^[™«l<5(i?-r[m(,)]) 

™{r) 

= e^^P[m(,)]5(i?-r[m(,)]) 
= e^^P[m(,)](5(i? + f[m(,)]) 



= e^P(-i?) . 

By integrating e~^P{R) = P{—R) over R, we get 

(e-'^[™(-)l) = 1 . (AlO) 

It follows from Jensen's inequality (e^) > e^^\ that 
(r[m(,)]) > 0. 

We now make a first choice of Wm,m'{Xt) in the back- 
ward dynamics ()A2p 

(All) 



W^"^ AXt) = W^"^ AXt) . 



In this case the backward dynamics is identical to the 
original one, except that the driving protocol is time re- 
versed. If we also choose the initial conditions of the 
backward dynamics to be the final conditions of the for- 
ward dynamics Pm(0) = Pm{T), using ()A6p and ([27|) . we 
find 



(A12) 
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The FT (j59j) previously derived using GFs now follows 
from Eq. (jAlOp . Using (|A9p , we also get the detailed 
form of the FT 



(A13) 



We now make a second choices of Wm,m'{^t) in (|A2p 
<L(At) = W^UX,)§^ . (A14) 



In the theory of M Es. (|A2|) with ()A14p is called the time 
reversal ME of ^ |29l, l3(|. We again choose the initial 
condition of the tilde dynamics to be the final conditions 
of the original dynamics Pm(0) = Pm{T). Using (|A7p 
with (|XT4|) and (|32|), we get 



(A15) 



The previously derived FT (|6Q|) follows now from (|AlQp . 
From Eq. ()A9p we find the detailed form of the FT 



P(A5na) 
P(-A5na) 



(A16) 



We can interpret the change in the total TEP during the 
to T time interval as the logarithm of the (forward) prob- 
ability that the driven system follows a given trajectory 
divided by the backward probability that the system, ini- 
tially in the final probability distribution of the forward 
evolution, and driven in a time reversed way compared 
to the forward evolution, follows the time-reversed tra- 
jectory. 

The nonadiabatic TEP is interpreted as the logarithm of 
the (forward) probability that the driven system follows 
a given trajectory divided by the backward probability 
that the system, initially in the final probability distri- 
bution of the forward evolution, and described by the 
time-reversed ME, follows the time-reversed trajectory. 
It should be noted that the backward ME (|A2p with 
(|A14p is different from the backward ME (|A2]) with 
(|Alip only for systems interacting with multiple reser- 
voirs which break the DBG. Only in this case ([59|) is 
different from (l6Ql). 
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FIG. 2: (Color online) (a) Driving protocol of the left lead 
chemical potential /ii(t) = /io + V{t) where V{t) follows 
(|97|) with Vi = -0.25, Vf = 0.5, c = 0.2, tm = 0.2 and 
fio — 0.5. The right lead chemical potential remains constant 
at fir = /io- (b) SoUd: The probability distribution of the 
dot obtained by solving the ME (|93p with initial condition 
po(0) = 0.4 and pi(0) = 0.6. Dotted: The adiabatic proba- 
bility distribution, (c) Decomposition of Stotif) according to 
([2T]) . (d) Decomposition of Sna(t) according to ([22)) . (e) De- 
composition of Snait) according to (|20p . (f) Decomposition 
of Stot{t) according to (|9]). 
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FIG. 3: (Color online) Five driving protocols for the left lead 
chemical potential /iz(t) = //o + y{t) where V{t) follows (|97p 
with Vi = -0.25, Vf = 0.5 and (ii) c = 0.2, tm = 75, (iii) 
c 0.05, tm 200, (iv) c 0.02, tm 300 (v) c 0.01, 
tm = 500. /ir(t) = /io = 0.5. (i) is the sudden switch limit 
l^i(t) = /io + ©(t)V/ [limit c ^ oo with tm = in (l97l) ]. 
The system is initially at steady state where po* =0.96 and 
pf = 0.04. In the final steady state pg* = 0.64 and pf = 0.36. 
a, b, c, d correspond to different measurement times. In all 
calculations /3 = 5, e = 1, = 0.2 and ar = 0.1. 
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FIG. 4: Probability distribution of the change in the system 
TEP for the protocols shown in Fig. [3] 
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FIG. 5: (Color online) Probability distributions of the change 
in the excess and driving TEP for the protocols in Fig. [3] All 
curves on the left column satisfy the FT (|60}. 
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FIG. 6: (Color online) Probability distributions of the change 
in the reservoir TEP for the driving protocols and measure- 
ment times shown in Fig. [3] All curves [except (i)a and (id) 
where the driving starts at the same time as the measurement] 
satisfy the FT {[6Z|. 




FIG. 7: (Color online) Probability distributions of the change 
in the adiabatic TEP for the driving protocols and measure- 
ment times shown in Fig. [3] All curves satisfy the FT (|61|) . 




FIG. 8: (Color online) Probability distributions of the change 
in the total TEP for the driving protocols and measurement 
times shown in Fig. O All curves satisfy the FT (|59)) . 
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FIG. 9: Summary of the two splittings of the total TEP in dif- 
ferent parts and of the type of transformations during which 
these parts are zero. SS-SS: transition between steady states. 
Eq-Eq: transition between equilibrium states. AdTrf: adia- 
batic transformation. TRSS: transient relaxation to steady 
state. TREq: transient relaxation to equilibrium. SS: steady 
state. Eq: equilibrium. <>= means that the ensemble 
average vanishes. 



